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Abstract: We derive a factorization theorem for Drell-Yan process at low qt using effec- 
tive field theory methods. In this theorem all the obtained quantities are gauge invariant 
and the special role of the soft function-and its subtraction thereof-is emphasized. We 
define transverse-momentum dependent parton distribution functions (TMDPDFs) which 
are free from light-cone singularities while all the Wilson lines are defined on-the-light-cone. 
We show explicitly to first order in ag that the partonic Feynman PDF can be obtained 
from the newly defined partonic TMDPDF by integrating over the transverse momentum 
of the parton inside the hadron. We obtain a resummed expression for the TMDPDF, 
and hence for the cross section, in impact parameter space. The universality of the newly 
defined matrix elements is established perturbatively to first order in Og ■ The factorization 
theorem is validated to first order in Os and also the gauge invariance between Feynman 
and light-cone gauges. 
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1 Introduction 

Physical observables with non-vanishing (or un-integrated) transverse-momentum depen- 
dence are specially important at hadron colliders. Those observables are relevant for the 
Higgs boson searches and also for proper interpretation of signals of physics "beyond the 
Standard Model". The interest in such observables goes back to the first decade imme- 
diately after establishing QCD as the fundamental theory of strong interactions [1-5]. 
Recently however there has been a much renewed interest in gy-differential cross sections 
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where hadrons are involved either in the initial states or in the final ones or in both [6-10]. 
The main issues of interest range from obtaining an appropriate factorization theorem and 
resumming large logarithmic corrections on one hand and to perform phenomenological 
predictions both for qq and gg production channels on the other. 

In this effort we re-examine the derivation of the factorization theorem for Drell-Yan 
(DY) heavy lepton pair production at small transverse momentum q^ and the proper 
definition of the non-perturbative matrix elements that arise in such factorization theorems. 
The region of interest is Aqcd <C ^ M, where M is the heavy lepton pair invariant mass. 
This topic was considered long time ago by Collins, Soper and Sterman [11] where the two 
notions of factorization and resummation of large logarithms (of the form a"ln"^(g|,/M^)) 
where systematically investigated. Those efforts yielded the well-known "CSS formalism". 
Here however we implement the techniques of soft-collinear effective theory (SCET) [12-15] 
to formally derive the factorization theorem for the gr-differential cross section. Within 
the framework of effective field theory, other efforts for g^-dependent observables were also 
considered in [6-9, 16]. 

Formal manipulations in SCET give us a factorized cross section for Drell-Yan at low 
qx, which pictorially speaking, looks as follows ^ 

= HiQ^flDPnifip) ® ® , (1.1) 

where = and fip is a factorization scale. Here H, F^^^^-^ and $ stand (respectively) 
for the hard part, the would-be (two) transverse-momentum dependent parton distribution 
functions (TMDPDFs) for the two collinear [n and n) directions of the incoming hadrons 
and the soft part. The above result might look familiar and in fact it resembles the one 
obtained by Ji, Ma and Yuan [17] for semi-inclusive deep-inelastic scattering (SIDIS) (with 
the relevant adjustments that need to be made when considering DY instead of SIDIS or 
vice versa). In Eq. (1.1) power corrections of the form {q^/Q"^)^ have been omitted. 

Explicit operator definitions for the various quantities in Eq. (1.1) will be given in 
the next Sections. However it is worthwhile at this stage to emphasize the important 
features implied in the derived factorization theorem, Eq. (1.1). The soft function ^{qr/Q) 
encodes the effects of emission of soft gluons into the final states with momenta that 
scale as Q{X,X,X) where A is a small parameter of order qr/Q- Those final state gluons 
(which hadronize with probability 1) are needed to kinematically balance the transverse 
momentum of the produced lepton pair. As we argue below, this function depends only on 
the transverse coordinates x± and the renormalization scale /i (which is implicitly assumed) . 
This feature of the soft function is consistent with the definition of the soft functions of Ji. 
et al. and Collins [17, 18]. The importance of such soft gluons was also acknowledged in [8] 
however due to the use of a special regulator, the "analytic regulator", their contribution 
vanishes in perturbation theory due to scaleless integrals. It is worth mentioning that 
in different regulators this is not the case and the soft contribution has to be included- 
explicitly-in the factorization theorem thus one obtains a regulator-independent theorem 
as the case should be. 

^The leptonic contribution is well-known and it is not shown in this section to simplify the notation. 
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In the effective tfieory approach soft and cohinear partons (with scahng of (5(1, A^, A) 
for n-colhnear or Q(A^, 1, A) for n-colhnear) are not allowed to interact simply because the 
collinear partons will be driven far off-shell. This is in contrast to ultra-soft and collinear 
interactions where ultra-soft gluons scale as Q{}? ,}? ,}?). However ultra-soft gluons are 
not relevant to the kinematical region of interest [6, 8] and will not be discussed further. 
The relevant framework to describe soft gluons (with A ~ qr/Q) interacting with collinear 
partons (with off-shellness Q^X'^ ^ ^qcd) named "SCET-gj^" in [6] and here we will 
adopt the same terminology. In SCET-gr the virtuality of the particles is of order q^, 
so it is different from SCET-II [19, 20], where the virtuality is of order Aq^Z)- SCET-II 
is needed in order to perform an operator product expansion (OPE) at the intermediate 
scale qT which would result in the appearance of the fully integrated PDFs. In both of 
these theories soft partons are decoupled from the collinear ones and their mere effect is 
manifested through the appearance of soft Wilson lines at the level of the matrix elements 
or Green's functions of the theory. 

Due to the fact that the soft function has a non- vanishing contribution in Eq. (1.1), 
then one needs immediately to consider the issue of double counting arising from the over- 
lapping regions of soft and collinear modes (when perturbative calculations are performed 
for the partonic versions of the hadronic matrix elements.) It turns out that this issue 
will dramatically affect the proper definition of the collinear matrix element (s), namely the 
TMDPDFs. In the traditional perturbative QCD framework the issue of double counting 
was treated through the notion of "soft subtraction" [21, 22]. In SCET, the analogous 
treatment was handled through "zero-bin subtractions" [20]. For sufficiently inclusive ob- 
servables (and at partonic threshold) an equivalence of the two notions was considered in 
[23-25]. Also in [6, 26] such equivalence was demonstrated up to C(a<j). In our case, we 
show in Sec. 9 that for certain IR regulators and in the kinematic region of interest, the 
equivalence of soft and zero-bin can also be established explicitly to first order in Og. 

Given the above and in order to cancel the overlapping contributions, the factorization 
theorem now reads 



where the small letter notation refers to the partonic versions of the collinear matrix ele- 
ments and to the partonic vacuum in the soft function (compared to the QCD one). The 
hatted symbols refer to the perturbative calculation of the collinear matrix elements that 
still include the contamination from the soft momentum region. Variations of the last 
result appeared in [17, 27] and also very recently in [18, 28] (see also [29], [10]). In SCET 
one could also consult [6, 23, 26, 30, 31]. 

Interestingly enough the last version of the factorization theorem, Eq. (1.2), is still 
problematic. Individually, the partonic collinear matrix elements and the soft function are 
plagued with un-regularized and un-canceled divergences which render them ill-defined. 
Those divergences show up perturbatively through integrals of the form: 



da = H{Qyfi^)[fn cl)-'] [fn 0"'] 



(1.2) 




(1.3) 
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which are manifestations of the so-called "light-cone singularities". Those divergences ap- 
pear, for certain IR regulators, also in the (standard) integrated PDFs, however they cancel 
when combining real and virtual contributions. This is not the case though for TMDPDF. 
Those light-cone divergences are a result of the fact that the Wilson lines (both soft and 
collinear) are defined along light-like trajectories thus allowing for gluons with infinite ra- 
pidities to be interacted with. To avoid such singularities, an old idea, due to Collins and 
Soper, is to tilt the Wilson lines thereby going off-the-light-cone. This trick was pursued 
in [17, 27] . More recently, Collins [28] argued that such regulator is necessary to sepa- 
rate ultraviolet (UV) and IR modes thus establishing two purposes: obtaining well-defined 
objects (free from un-regularized divergences) and a complete factorization of momentum 
modes. In the light of the above arguments, one needs to define a new set of collinear 
matrix elements as follows 

fntjnt) . 
Jnt{nt) f-T- 1 

where the subscript t stands for "tilted" Wilson lines which are no longer light-like. This 
is true for collinear and soft ones as well. With this, the factorization theorem takes the 
form 

da = H{Q'^/f?)jn,®jnt . (1.5) 

In this work we take a different path. We show that all IR divergences, namely the 
soft and collinear ones appearing in a massless gauge theory, as well as the light-cone 
singularities can still be regularized while keeping all the Wilson lines defined on-the-light- 
cone. When going off-the-light-cone one introduces the (^-parameter: C, = (Prit)'^ / (rit)^ 
where P stands for the incoming hadron momentum. This parameter complicates the 
phenomenological studies since, among other things, it will affect the evolution of the 
hadronic matrix elements. However when staying on-the-light-cone, the evolution of the 
TMDPDF will be governed only by the factorization scale ii. Second-and on the technical 
side-the non-vanishing small components of rit and nt introduce small contributions (in 
powers of effective theory parameter A) that violate the power-counting of that theory 
unless some ad-hoc relations are imposed between the small and large components of the 
tilted vectors. It is also not so clear how one can relate the TMDPDF with the integrated 
PDF when going off-the-light-cone. Moreover, staying on the light-cone is much more 
compelling when one considers computing, say, the TMDPDF and its anomalous dimension 
in light-cone gauge. When choosing this gauge then going off-the-light-cone is completely 
awkward. Those considerations motivate us to stay on-the-light-cone. When doing so one 
gets 

da = H{Q^/fl^)jn®jn. (1.6) 

The above result is still an intermediate step towards getting the final factorization 
theorem. However, we will define our TMDPDFs based on it, where j„(n) will be the par- 
tonic TMDPDFs, from which one can easily get the hadronic ones. An extended discussion 
of the "on-the-light-cone TMDPDF" and its properties will be given in the Sections below 
and it forms the basic contribution of this effort. 



-4- 



Given that the TMDPDFs, 

jn{n)^ include soft contribution then they become dependent 
on the perturbative intermediate scale qt-, thus a further step of factorization is needed. 
This is achieved by performing an operator product expansion (OPE) in impact parameter 
space in the region b <^ ^qcd^ where the TMDPDFs are matched onto the PDFs Qn(n)- 
Once performing the OPE in impact parameter space we get ^ 



dT, = H{Q'^/fi'^) Cn{x;b,Q,fj.) (g) Qn{x;fi) Cniz;b,Q, fi) ® Qn{z; fi) 



(1.7) 



Notice that Cn{n) still have an explicit Q^-dependence. This dependence is harmless in 
the sense of factorization of Lorentz invariant scales, since H and Cn{n) ^-re both pertur- 
bative while Qn{h) ^'I's non-perturbative. However this dependence asks for resummation 
of logarithms of Q'^ j ^ once [i is chosen to be much smaller that Q. The extraction of 
Q^-dependence of Cn{n) a^id its resummation thereof is discussed in Sec. 5 and the final 
result for the cross section is 



4e-27B 
4e-27B 



(1.8) 



This result allows the resummation also for large logarithms of Aqcob to be performed by 
simple running between different scales. 

Another novel feature of our derived factorization theorem is gauge invariance. Re- 
cently it was shown [32, 33] that SCET, as was traditionally formulated, has to be adjusted 
by the inclusion of transverse Wilson lines, Ts, so as to render the basic building blocks 
and the Lagrangian of SCET gauge invariant under regular and singular gauges. This has 
the powerful result that all the derived physical quantities (appearing for example in the 
factorization theorem Eq. (1.1) or the likes) are gauge invariant and no transverse gauge 
links need to be invoked by hand in the aftermath. This derivation allows one to consider, 
for example, the subtracted TMDPDF in covariant gauge, say Feynman gauge, and a sin- 
gular gauge, say light-cone gauge. We have carried out such computations and found, as 
expected, full agreement to hold at first order in the strong coupling a^. 

It is clear that the soft function $ (or </>) connects two different collinear sectors. 
This might spoil the single-collinearity notion inherited in the standard definitions of the 
PDFs or TMDPDFs thus the universal features of such quantities might be jeopardized. 
We discuss this issue extensively and establish the universality of the TMDPDF given in 
Eq. (1.4), in agreement with [28]. We also establish, for the first time-and to first order in 
Og-that when integrating over transverse momentum, the integrated PDF can be recovered 
from the TMDPDF. 

This paper is organized as follows. In Section 2 we establish the factorization theo- 
rem for DY at small qt and the proper definition of the TMDPDF on-the-light-cone. In 



^Notice that the convolution in Eq. (1.6) is in momentum space with respect to k„± and kfij_ while the 
convolution in Eq. (1.7) is in the Bjorken variables x and z. 
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Section 3 we calculate jn{n) its anomalous dimension to first order in a^, and discuss 
the absence of rapidity logarithms in the TMDPDF. In Section 4 we show explicitly how 
to obtain the standard PDF from the TMDPDF to first order in by integrating over 
the transverse momentum of the parton. In Section 5 we discuss the Q^-dependence of the 
intermediate matching coefficients and thus establish the complete factorization of scales. 
We also show the resummed TMDPDF (and hence the cross section), give its anomalous 
dimension at second order in a^. In Section 6 we demonstrate that the TMDPDF attains 
the same order as result both in DY and DIS kinematics. In Section 7 we establish the 
validity of the factorization theorem to first order in as ■ The gauge invariance of our result 
is discussed in Section 8 and Section 9 shows the equivalence of zero-bin and soft function 
subtraction. Then we conclude in Section 10. In Appendix A we show the calculation of 
the matching of the TMDPDF onto the PDF with 5-regulator for DIS and DY kinematics. 
In Appendix B we present the calculation of the Quark Form Factor in full QCD with 
(5-regulator for DIS and DY kinematics. 

2 The Factorization of Drell-Yan at Small qx 

Let the momenta of the two incoming partons initiating the hard reaction be p and p. A 
general vector v is decomposed as: = h- + n ■ +v± = v~^^ +v~^ + v±, where 
n = (1,0,0,1), h = (1,0,0,-1). ^ We denote v = \v_l\, and particularly qx = \q±\- The 
momentum scalings of the n-collinear and n-collinear were given in the previous Section. 

Together these modes give the momentum scaling of the outgoing photon: Q{1, 1, A). 

The initial form of the cross section is 

'^^^W'sWf^ ^ / Ae-*^-^(-5M.)(^^i(^',^i)^^2(P,a2)|J'^t(y)j-(0)|iV,(p,^^)iV2(P,a2)), 
J'^ = ^e,^7^V', (2.1) 

where J is the electromagnetic current and Cg is the quark electric charge. P and P 
correspond to the hadrons momenta and s = {P + P)"^. Note that the scaling of the 
position variable y in Eq. (2.1) is y ~ 1/(5(1, 1, 1/A) and we will make use of this below. 
The full QCD current is then matched onto the SCET-g^ one 

r = C(QV^2) ^ e.XnSVl^Slxn , (2.2) 
<? 

where in SCET the n-collinear and n-collinear (or anticollinear) fields are described by 

^ Notice that our convention for the light-cone coordinates is opposite to the one used in the standard 
SCET Kterature. 
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Xn(n) = W^J(ii)Cn(n)- For DY kinematics we have 



n{n) 



Wn{x) = Pexp 
Tnix) = Pexp 
Tn{x) = Pexp 



I'd 



ds n ■ An{x + sn) 



dr l± ■ An±{x~^ ,00 ,x± + 1±t) 

1 

dr l± • An_L{oo~^ , x~ ,x± + /_|_r) 



(2.3) 



Wfi can be obtained from Wn by n -H- n and P ■v^ P, where P{P) stands for path (anti- 
path) ordering. The transverse Wilson hnes are essential to insure gauge invariance of 
Xn[n) among regular and singular gauges [33]. 

The soft Wilson lines and their associated transverse Wilson lines are given by 



'^n(n) 



^ sn(sn) *-'n(n) ) 



Sn{x) = Pexp 
Tsn{x) = Pexp 
Tsfiix) = Pexp 



ig 
ig 



ds n ■ As(x + sn) 



00 


-00 




dr l± • As±{oo~^ ,0 ,x± + 1±t) 
dr l± • As±{0~^ , oo~ ,x± + 1±t) 



(2.4) 



where Tg^(^sf^^ appears for the gauge choice n • yl^ = (n • = 0) and Sn can be obtained 
from 5„ by n -H- n and P P. In the SCET literature there are different ways for 
obtaining the appropriate soft and collinear Wilson lines. However, one can also start from 
the full QCD vertex diagram and then take the soft or the collinear limit of the virtual 
gluon loop momentum. The resulting vertices obtained can unambiguously determine the 
soft and collinear Wilson lines in the effective theory. The above definitions of the Wilson 
lines are compatible with the QCD soft and collinear limits for time-like (DY) virtualities. 
In Sec. 9 we present the Wilson lines relevant for space-like (DIS) kinematics and their 
derivation follows the same argument of taking the soft and collinear limits of QCD. 

Using Fierz transformations and averaging over nucleons spins, the hadronic matrix 
element in Eq. (2.1) can be casted in the form 



- ( A^l (P, ai ) iV2 (P, (72 ) I J'^T (y ) (0) I AT, (p, ) (P , (72 ) ) 



|C(QV^)|2^e2 — (iVi(P,c7i)iV2(P,^2)| (x^(y)ix*i(0)) (xn(y)^Xn(O) 



X Tr 



T(5t(y)5^(y))T(4(0)5„(0))l |iVi(P, c7i)iV2(P, (72)). 



(2.5) 



Since the n-collinear, n-collinear and soft fields act on different Hilbert spaces one 
can disentangle the Hilbert space itself into a direct product of three distinct Hilbert 
spaces [12, 25]. The collinear, anticollinear and soft fields obey different Lagrangians 
which are opportunely multipole expanded [12], however the multipole expansion of these 
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Lagrangians does not affect the "^''-dependence of the fields in Eq. (2.5) (because there 
are no interactions among soft and colhnear fields). Due to these arguments, one can then 
write the cross section as 

dS = ^^^ Jd^e-"^-y H{Qyf^')Y.elFr.{y)F^{y)^y), (2.6) 
where HiQ^fi^) = |C7(Q2/^^)P and 

Fn{y) = \ Y.^N,{P, CTi)| Xn{y) I Xn(0) |iVl(P, di)) , 
Fn{y) = \ E(^2(P, 02)\ Xm \ Xn{y) |A^2(P, ^2)) , 

^{y) = k^\^T{Sl^Sl\{y)T{sVsl\m^) ■ (2-7) 

We now Taylor expand Eq. (2.6) in the physical limit that we are interested in. The 
photon is hard with momentum q Qil, 1, A), so the in exponent e"*''^ in Eq. (2.6) one 
has y -^{1, 1, 1/A) as mentioned before. On the other hand the scaling of the derivatives 
of the n-collinear, anticollinear and soft terms are clearly the same as their respective 
momentum scalings. Combining this with the scaling of y, the leading term (0(1)) of the 
cross section reads 



2 
Q 

q 



SNcQ^s (27r)4 

xF„(0+,2/-,yl)i^n(y+,0-,y^)$(0+,0-,yl) +0{\). 



It should be immediately noted that if one had considered ultra-soft scaling Q{X'^, A^, A^) 
instead of the soft one in then after the Taylor expansion <I> would be exactly 1 to all 
orders in perturbation theory. This is the case that was considered in [8]. The fact that 
the soft function, $, depends only on the transverse coordinates is of crucial importance. 

In the rest of this Section we will consider the leading order contribution to the partonic 
version of the cross section Eq. (2.8), namely 

3Ncq^ 2(27r)4 w ; . 



1 



with 



j d'^kn±d'^kn±d'^ks± 6^^\q± - kn± - kn_L - ksi_) fn{x; kn±) fn{z; kni_) (f>{ks±) ■ (2i 



2 J (27r)= 

m±) = J ^e-"-'="'^(0+,0-,f^) . (2. 
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Discussion of the TMDPDF: 

The problem with the above factorization formula is that in the 0{as) calculation of 
fn{h) ^-iid (j) there are still mixed UV and IR divergences which complicate both the renor- 
malization procedure and the non-perturbative interpretation of such quantities. However, 
when considering the following combinations 



jn{z] kn±) 



1 f dr (Pf± 



1 f dr^(fr 



-i{\r xp+-r^-kn±) f (n+ 



/n(O+,r-,fx)V0(O+,O-,f^) 



(2.10) 



it turns out that those quantities are free from such mixed divergences. This is shown 
explicitly to hold to 0{as) in the next Section. jn{n) in the last equation are our definition 
of the TMDPDFs. 

As mentioned in the Introduction, the equivalence of the zero-bin and soft subtractions 
is a subtle issue, since it depends on the kinematics and the IR regulator that is used. In 
Sec 9 we show this equivalence to first order in a^, which leads us to the following relation 
between / and /, 



/n(0+,r",r-l) 



A(0+,r-,fl) 
(/.(0+,0-,r-l) 



where 



/n(r+,0-,f=l) 



/„(0+,r-,r-l) = {p\ [^nW^\ (0+,y",yx)| [W^^ (0) |p)|,b included : 



/n(r+,0-,fl) 



(2.11) 



/n(r+,0-,f=l) = { 
</.(0+,0-,f=l) = (0|Tr 



CnWf ] (0) I [Wr^ {y^,0~,y±) \P)U included : 



(0+,0-,y^) 



qT] qT 



(0)|0) . 



(2.12) 



Then, we can write the TMDPDFs in the following way, 



1 f dr d'^r± 



A(0+,r-,rl) 



Jn{Z 



(27r)3 ^ v/'A(0+,0-,r-l) 

1 r rfr+dVl ^^(^,+,p--,., .fc, , ) /n(r+,0",r±) 
'^"^^-27 (2vr)3 v/</'(0+,0",r-l) 



(2.13) 



where the square root of the soft function is subtracted from the naive collinear matrix 
element. 

Thus it is compelling to re-cast the factorization theorem in the following form 
Atto^ dxdzd^qj_ 



-H{Q'/f,')J2e 



da 



which is formally in agreement with the line of argument of Refs. [18, 28]. The TMD- 
PDFs 

jn{h) defined in general in Eq. (2.10), but in the following we take the result in 



SNcq"^ 2(27r)4 

X I d^knLd^kn±5'''^\q^-kn±-kni_)in{x]K (2.14) 
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Eq. (2.13), which apphes for our particular kinematical regime (perturbative qt and away 
from threshold) and for the set of IR regulators implemented below. 

The last result is still not the final form of the factorization theorem and the TMDPDFs 
still have to be refactorized. In the effective theory approach this corresponds to a second 
step matching of SCET-g^ that describes the physics at the intermediate scale qt ^ ^qcd 
with SCET-II that captures the non-perturbative physics at the hadronic scale Aqcd- 

The refactorization of the TMDPDF is essential since the naive collinear fn{h) and 
the soft contribution (j) that enter in the definition of jn{h) live at the intermediate scale 
qt, consistent with their construction in SCET-qt. Since qt is perturbative, its conjugate 
coordinate, the impact parameter b, is small enough to perform an OPE in the impact 
parameter space. Moreover in this space the IR structure becomes manifest with the 
appearance of IR poles in dimensional regularization. Obviously, the first term in the OPE 
would be just the standard Feynman PDF, and the Wilson coefficient would be the term 
that sums all the large logs between Aqcd and qx (see Refs. [4, 5, 11] and more recently 
using SCET Refs. [6, 8, 34, 42]). Then, given the following OPE (and an analogous for n) 

jn{x; b^,fi) = / —Cn -; b^, fi Qn{x'; fi) + Oib^Alco) , (2.15) 

Jx X \X J 

where 

~jn{x;b^,ll)= [ d^kn±e'''"^-''^jn{x;knL,fJ-), (2.16) 



the factorization theorem takes the form 



47ra^ dxdzd'^q± sr^ 2 f d?b^ ~iqj_-b, dz' 



3iV,(?2 2(27r)4 ^ (27r)2 

x//(QV/"')C'n(^;^±,Q,/i) Cn(^^-A,Q,f?j Qn{x';f^)Qn{z';fi). (2.17) 

In this effort and for simplicity of presentation we will not consider the contribution 
coming from a gluon splitting into two quarks. This contribution is certainly vital for the 
final result of the DY cross-section. Here however we are mainly interested in studying the 
TMDPDF of a quark in a quark. Henceforth we will refer to this quantity simply as the 
"TMDPDF" and it can be easily checked that all the results below are not affected by this 
omission. 

The above result is one of the main results of this paper and it holds to all orders in 
perturbation theory. It is worthy to notice the separation of scales: the hard matching 
coefficient lives at scale Q, the matching coefficients at the intermediate scale live at 1/5 ~ 
qt, and finally the PDFs live at the hadronic scale Aqcd- 

As we show below in Sec. 5, C'n(n) have a subtle Q^-dependence which at first sight 
might spoil the scale factorization, however this dependence can be extracted and expo- 
nentiated thus putting it under control. 
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3 The TMDPDF On-The-Light-Cone 



In this section we compute the TMDPDF in Eq. (2.13) to first order in Os while regular- 
izing the UV divergences in dimensional regularization (DR) and in MS scheme. The IR 
divergences as well as the light-cone ones will be regularized by the 5-regulator introduced 
below. We also present a calculation of the virtual part of the TMDPDF in pure DR. 
This is shown in Subsection 3.5. Notice that individual Feynman diagrams with real gluon 
emission cannot be obtained in pure DR due to light-cone singularities and an additional 
regulator is needed. When such diagrams are added up, the logarithmic dependence on this 
regulator cancels out which means, at least to first order in a^, that real gluon contribution 
is free from light-cone divergences. More discussion on this is given below. 

We write the poles of the fermion propagators with a real and positive parameters 



(p + A;)2 + iO {p + k)^ + iA-' {p + k)'^ + iO {p + ky + iA+ 



(3.1) 



The above prescription applies as well to the fermion propagators in SCET. The cor- 
responding pole-shifting for collinear and soft Wilson lines goes as follows 



1 



1 



1 



1 



k+ ±iO k+± i6+ ' k- ±iO fc" ± i6- ' 
where 5^ are related with through the large components of the collinear fields 



A' 



(3.2) 



(3.3) 



The 5-regulator resembles the one used in [26]. Notice that the essential feature of the 
(5-regulator is the coherency between the regulators of fermion propagators and the corre- 
sponding Wilson lines. 

Expanding Eq. (2.13) to first order in Os one finds 



jn{x', kn_L, Q, /i) 



1 fdC d'^^^ ^-i{U-xp+-i^-K^) 



2 J (2^)2 



-e ^z'' 



/nO + ( /nl — ^/nO*^! 



0{a 



si ; 



(3.4) 



where the numerical subscripts denote the order in the Ug expansion. The collinear matrix 
element at tree level is 



/no = (Pie„(r ,o+,e±) f um = e*^^^«>+ , 



(3.5) 



so the expansion of j„ up to order is 



jn{x;kn±,Q,n) = S{1 



+ 



x)(5(2)(fe 



(27r)3 



J'l 



nl 



x) 



(2vr)2 



(3.6) 
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(a) 



(b) 



Figure 1. Virtual corrections for the collinear matrix element. The black blobs represent the 
collinear Wilson lines W in Feynman gauge or the T Wilson lines in light-cone gauge. Curly 
propagators with a line stand for collinear gluons. "h.c." stands for Hermitian conjugate. 
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+ h.c. 



(a) 
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(c) 



Figure 2. Virtual corrections for the soft function. Double lines represent the soft Wilson lines, 
S. "h.c."' stands for Hermitian conjugate. 

3.1 Virtual Diagrams 

The diagrams in figs. (1) and (2) give collinear and soft virtual contributions respectively to 
jn- The Wave Function Renormalization (WFR) diagram (la) and its Hermitian conjugate 
give 



1 1, 1 

h -In— 4 h - 

2euv 2 -iA- 4 



+ h.c. 



(3.7) 



The W Wilson line tadpole diagram, (lb), is identically 0, since = 0. Diagram (Ic) and 
its Hermitian conjugate give 



h.c. 



2tt 



6{1 - x)5^^Hkn±) 



2,5+ 2 
-In— + 



,5+A' 



A- 



euv 



In^— ^ - 21n^ + In^— + 2 - — 

/i^ 12 

(3.8) 



The contribution of diagrams (2a) and (2b) is zero, since (2a) is proportional to = 
and (2b) to = 0. The dia gram (2c) and its Hermitian conjugate give 



2e 



1 



{27r)d [k+ - i5+] [k~ + iS-] [A:2 + iO] 



+ h.c. 
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,5+6- 



— + Y 



(3.9) 
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Figure 3. Real gluon contributions for fn{n) 



The virtual part of the TMDPDF at 0{as) using the relation in Eq. (3.3) is 



1 1 

^ + 



'UV 



3, A 1, 
-m- 



2 //^ 



-In^ 



A2 



+ W 



+ 



7 -k' 



3 



(3.10) 



where we have set A"*" = A = A. Analogously the anticollinear one is 



2-K 

2 /i2 2 



1 1 

^ + 



'UV 



A2 



.A 7 7r2 
/i^ 4 3 



(3.11) 



As mentioned earlier, individual contributions to j^i{ni) have mixed UV and IR divergences, 
however j'" itself is free from them. 

3.2 Real Diagrams 

The relevant diag icims for tli6 rG3,l part of jVi are shown in figs. (3) and (4). The diagram (3a) 
gives 



X 5 ((1 - x)p+ - k+) 5^^\k± + kn±) 



-^^{l - e){l - x) — 



(2vr)^ 



2 ' 



(3.12) 



where the e-dependence and the deltas are kept for later convenience. The e will be needed 
when transforming to impact parameter space in pure DR, and hence putting all deltas to 
zero. The deltas will be needed when going to impact parameter space in 4 dimensions, 
using them to regulate k± — )■ 0. The sum of diagram (3b) and its Hermitian conjugate (3c) 
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Figure 4. Real gluon contributions for the Soft function. 
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?(3fe+3c) 
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+ -k+ 



^ ^ J {2Try ^ ' ^ ' [k+ + i6+][{p - kY + 
x5({l- x)p+ - A;+) (5(2) (A?^ + 4^) + h.c. 



2asCF_ 
(27r)2-2^ 



(1 - x) +i6+/p^ 



iA-{l-x)_ 



+ /i.e. , 



(3.13) 



Diagram (3d) is zero, since it is proportional to = 0. 

For the real emission of soft gluons, the diagrams (4a) and (4d) are zero, since they 
are proportional to = and = respectively. The diagram (4b) and its Hermitian 
conjugate (4c) give 



,(46+4c) 



i-jdt. 1 

1 , 5+5- 
-m- 



(2vr) 



2-2e 



kn±?-5+5- 



\k , |2 



h.c. 
(3.14) 



To get the real part of the TMDPDF we need to put together all the contributions above 
into Eq. (3.6). In order to do that we set the deltas to zero unless they are necessary to 
regulate any divergence, i.e., we keep them in the logs but not for the limit \k±\ — ^ 0. The 
result, setting = A, is 



2asCF 1 



[l-e){l-x) + 



+ o(l — xjm 



(1-x) 



k 1 |2 



(3.15) 



where the InA has cancelled in the combination of Eqs. (3.13) and (3.14), using the rela- 
tion (A. 5) to express the first one in terms of distributions. A similar result is obtained for 
the anticollinear TMDPDF 



Jul 



2asCF 1 

W) 



2-2e 



h- 1 |2 



{l-e){l-z) + 



+ o(l — 2;)ln 



(1-z) 



h- 1 |2 



(3.16) 



where the superscript r stands for "real" gluon contributions. Our results for the real gluon 
emission agree with the ones in [9, 17]. 
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3.3 Where Are the Rapidity Divergencies? 

As can be seen from Eqs. (3.15-3.16) the contribution to the TMDPDF from real gluon 
emission does not include any logarithmic dependence on 5's after a fine cancelation of those 
logarithms between the naive collinear and the square root of the soft contribution. Clearly 
however this regulator is needed to calculate individual Feynman diagrams in order to avoid 
light-cone singularities. Now let us consider the contribution from virtual diagrams given in 
Eqs. (3.10-3.11). In obtaining those results there is also a fine cancelation of the pieces with 
mixed divergences appearing in both terms of the soft and naive collinear contributions. 
However the remaining A-dependence can be safely interpreted as pure IR one, namely 
soft or collinear, and not as rapidity divergence. The reason is simple. In Appendix B 
we have calculated the quark form factor in full QCD using the same 5-regulator. Since 
QCD has no rapidity divergences (or in other words there are no Wilson lines) then all the 
A-dependence is pure IR. Those IR terms (see Eq. (B.5)) are exactly the ones appearing in 
Eqs. (3.10-3.11). Among other things this fact allows for QCD to be matched onto SCET 
at the higher scale Q thus obtaining a matching coefficient, i?, which is free from any IR 
dependence. Moreover the logarithms of A's in the virtual contribution need not to be 
resummed. So we conclude that our TMDPDF is free from any rapidity divergences. 

To strengthen the above statement we mention that one can also obtain a A-free 
virtual contributions by simply calculating them in pure DR, as we show below. Also 
in that case there is a fine cancelation of mixed poles (UV and IR) between the naive 
collinear and the square root of the soft function and the virtual contributions contain 
both UV poles and IR ones. The IR ones will again exactly match the ones of QCD for the 
quark form factor. Then the complete first order result for the TMDPDF can be obtained 
without any 5-dependence. In other words it is free from any rapidity divergences. It is 
interesting to notice that also when using the regulator introduced in [34] (the z^-regulator) 
our observation above is still valid. In [34] the zero-bin contribution is zero to all orders 
in perturbation theory. However the soft function defined in that work, which at the 
operator level is identical to ours, is non-vanishing and z^-dependent. In this case we resort 
to our defining formula for the TMDPDF, namely Eq. (2.10). Taking the results from 
[34] (Eqs. (5.50, 5.51) for collinear and Eq. (5.62) for the soft contribution) one can easily 
verify that the TMDPDF, calculated with z/-regulator is also independent of the rapidity 
divergence regulator v. Although the above discussion relies on first order calculations in 
a^) it holds to all orders in perturbation theory as we explain below. 

The full QCD calculation of the partonic cross section can be performed while staying 
on-the-light-cone and using pure DR to regularize both the UV and IR divergences. When 
matching QCD and SCET at the higher scale Q the hard matching coefficient should 
be independent of any IR regulator to all orders in perturbation theory (as long as the 
assumption that SCET captures all IR behavior of QCD is still intact).^ Given the relation 
in Eq. (2.14), it is clear that jn(n) should also be free from any rapidity divergence regulator 

In [17, 27] the QCD matrix element was calculated on-the-light-cone, while the SCET ones where 
calculated off-the-light-cone. Thus, the hard matching coefhcient is dependent on the ofl-the-light-coneness 
parameter p. 
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since the full QCD result is. 

Consider the virtual contribution to the TMDPDF given in Eq. (3.10) calculated with 
the 5-regulator. Since the real gluon contribution (see Eq. (3.15)) is free from that regulator, 
then the TMDPDF does depend on it and this dependence persists also at higher orders 
in perturbation theory. However the important thing to notice is that the anomalous 
dimension 7„ of the TMDPDF and its matching coefficient C„ onto the PDF, which are 
the only perturbatively calculable physical quantities relevant to the hadronic TMDPDF 
(in the kinematical regime we are interested in), are still 5-fiee. This is due to the fact 
that in Eq. (3.10) we managed to separate the UV from all the IR divergences including 
the light-cone ones. In Section 4 and in Appendix A we show this explicitly to first order 
in as- The derivation of the anomalous dimension of the TMDPDF is given in Sees. 3.4 
and 3.5. 

As mentioned in the Introduction, one of the prominent methods to regularize light- 
cone singularities is to go off-the-light-cone where we tilt the light-cone vectors: (n, n) — 
(f,u), where ^ 0, v'^ ^ 0. The diagrams with real gluon emission that have light- 
cone singularities and need to be regularized are (3b, 3c, 4b, 4c). From the naive collinear 
contribution we get 



?(36+3c) _ 2asCF 1 



2^ M 
+ o(l — xjln- 



and from the soft contribution we have 

,(46+4c) _ 2asCF 1 



(3.17) 



- V , (3.18) 



where p = {v~ /v'^) = {v'^ /v^) and goes to infinity in the light-cone limit. Eqs. (3.17- 
3.18) are consistent with the results found in [17]. Combining these two results according 
to Eq. (3.15) with "J ^ given in Eq. (3.12), we get again the result in the second line 
of Eq. (3.15), where all the /9-dependence cancels in exactly the same manner as the 5- 
dependence was canceled. The last observation differs from the one found by Collins [28]. 
The reason is that in Collins' work, the contribution from the collinear matrix element 
to the TMDPDF is treated while staying on-the-light-cone and only the soft function is 
"tilted" (see Eq. (10.136) in [28]). In [17] (see also [27]) ah matrix elements, soft and 
collinear, are taken off-the-light-cone in a consistent manner. In this case, as we have 
shown, the p-dependence cancels in the TMDPDF. 

The authors in [34] have a different definition of the TMDPDF than ours in the sense 
that the soft function is not included in that definition. Moreover, since they identify 
the collinear matrix element as the TMDPDF, it contains rapidity divergencies, which are 
interpreted as UV ones (because of the regulator they use), and are then resummed using 
usual RG techniques. 

3.4 Anomalous Dimension of the TMDPDF jn{n) 

In order to compute the anomalous dimension of the TMDPDF to 0{as) one needs to 
consider only its virtual contributions since, as we have seen in the previous sub-section. 
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all real contributions are UV-finite. Prom Eq. (3.10) the counterterm for the n-collinear 
TMDPDF is 



1 



2-K 



1 1 



'UV 



(3.19) 



and the corresponding anomalous dimension is 
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1 dZn 1 dZn . 2^^ 

+ —- — (-2eas + 0(aj) 



din// Zn dlnn Z^ dag 



27r 



3 + 21n^ 



(3.20) 



For the n-collinear sector we have, analogously, Zn, from which we get 



7nl 



dlnZn 
dln/i 



<XsCf_ 
27r 



3 + 21n 



(3.21) 



Now we compare our result for the anomalous dimension with the one of Collins [28] 



^nl,Collins 



2tt 



3 + 21n 



2(p^ 



(3.22) 



where y„ is a measure of "off-light-coneness" which sets a lower bound on the rapidity of 
gluons. Collins uses tilted vectors along the trajectories of the incoming hadrons nt and hf 
which are parameterized as nt = (1, — e~^^", 0^) and ht = (— e^^", 1, Oj_). Thus it is clear 
that in order, for example, to recover the light-cone limit one has to consider two different 
limits: ?/„ — >■ oo for rit and y„ — )• — oo for ht- This means that the light-cone limit in 
Ini.Coiiins is Unattainable. Moreover, in Ref. [28], the ^-dependence cannot be eliminated 
by any single natural choice. In our case, however, we see that by setting one regulator, 
A = A"*" = A~, the anomalous dimension becomes independent of that regulator and it is 
free from rapidity divergences. 

Moreover, the evolution of the TMDPDF can be governed only by the //^-evolution 
and no additional parameters are needed. 

In Sec. 5.1 we give the AD of the TMDPDF at second order in a^. 

3.5 Anomalous Dimension On-The-Light-Cone in Dimensional Regularization 

It is possible to calculate the anomalous dimension of the TMDPDF using pure DR as 
in Ref. [23]. At one- loop, we need only to consider the virtual contributions given in 
diagrams (la), (Ic) and (2c). All the rest vanish identically due to light-like Wilson lines. 
For diagram (la) (without its Hermitian conjugate) we have 



4vr Veuv eii 



em) 



(3.23) 



^The factor of 2 difference between Eqs. (3.20, 3.21) and Eq. (3.22) in the logarithmic term is due to 
different conventions for the light-cone components. 
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and for diagram (Ic), 



fnl 



(Ic) _ asCp 



An 



6{1 - X)5''^\kr, 
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eiR 



+ 



2 



(3.24) 

Notice that in this regularization scheme the energy scale inside the logs is fixed noting that 
is the only relevant scale in the virtual part of the TMDPDF. Thus the scale inside the 
logs is equal to — where k = Q'^/{p'^)'^ and it is required to remove the dimensional 
ambiguity in integrals of the form: /q°° dt The soft function, diagram (2c), gives 



,,(2c) 



2tt 





■ 1 


1 " 










.euv 


eiR. 



(3.25) 
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Taking into account the Hermitian conjugate diagrams, the total virtual contribution to 
the TMDPDF is 



1 



21n 
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2euv 



1 21n 

— + — 
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(3.26) 

Prom the result for j^^ one can easily identify the counter-term Zn needed to cancel the 
UV divergences. Defining 7^ = one gets 



7nl 



27r 



3 + 21n 



27r 



3 + 21n 



(3.27) 



which agrees with Eq. (3.20). 

In Eq. (3.26) we again notice that there are no mixed UV and IR divergences as was 
observed in Eqs. (3.10, 3.11). It should be noted that if one had subtracted the complete 
soft function from the collinear part (and not the square root of it) then there would be 
mixed UV and IR poles and those mixed poles would not cancel even after including the 
contribution from real gluon emission. This would definitely prevent such quantity from 
being an acceptable definition of TMDPDF. 



4 From TMDPDF to Integrated PDF 

We recall from Section 2 the OPE of the TMDPDF onto the PDF 

jn{x;b±,Q,f^) -- 



(4.1) 



where 



Q^{x;fi) = ^ I ^e-^t?' -P+(p|^„(0+,y-,0x)|xU0+,0-,0x)|p)|zb included. (4.2) 



1 / ^e-^h 

In this section we compute Cn to first order in Og and also establish to the same order the 
following relation in d = 2 — 2e 

/i^^ / d'^kn±jn{h){x;kn±,Q,fJ,) = Qn{x; fi) , (4.3) 
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which is expected to hold to all orders in perturbation theory for bare quantities. We 
comment more on this below. 

In impact parameter space the virtual and real parts of the TMDPDF are obtained 
from Eqs. (3.26, 3.15) by Fourier transformation. The results are 
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where Lt = ln(/i^6^e^'^^ /4) and 



(4.4) 
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(1-x). 



+ (l-x) + -<5(l 



(4.5) 



which is the one-loop quark splitting function of a quark in a quark. We have used pure 
DR and dropped the A~ in the real diagrams. In Appendix A this calculation is done 
while keeping the A~ to regulate k^x — ?■ when going to the impact parameter space. 

The complete renormalized TMDPDF in impact parameter space, to first order in a^, 
becomes 



jn{x] b±, Q, /i) = 5{1 - x) + 
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Given that to first order in Og the renormalized PDF is 

asCp 



Qn{x; fx) = 5{1 - x) + 
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(4.7) 



we extract from Eq. (4.1) the matching coefficient at first order in Og, 

OisCp 



Cn{x; b, Q, fx) = 5(1 - x) + 
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-Lt + In^Lr + 



TT 

12 



(4.8) 



At this stage it is worth noticing the appearance of ln(Q^///^) at the matching coefficient. 
From the above result, we can see that by a proper choice of the scale fi = fii = (2e~'^^/6), 
we eliminate this logarithm since LT{fxi) = 0. However at this order in perturbation theory 
this cancelation is accidental and it does not persist at higher orders. In subsection 5 we 
discuss the appearance of ln(Q^/^^) at an arbitrary order in perturbation theory and how 
to handle them. 

It has been a matter of debate whether the PDF can be obtained from the naive 
collinear matrix element, fn(n)i by simple integration over the transverse momentum k± 
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conjugate to b± (see, e.g., [17, 28, 45] and more recently [29, 46]). In principle, and as a 
consistency check of the partonic definitions of such quantities, this should be the case; 
however, and interestingly enough, this was never established before. Below we establish 
that even with the inclusion of the soft function in the definition of the TMDPDF we are 
able to recover the PDF in a straightforward manner. It should be mentioned that going 
from TMDPDF to PDF can only be obtained when considering bare quantities. In other 
words if we consider the renormalized TMDPDF and then integrate over the transverse 
momentum we cannot recover the renormalized PDF since the integration over the trans- 
verse momentum introduces new UV divergences that need to be renormalized in their 
own turn. Stated differently, it is the lack of interchangeability between integration and 
subtraction (of UV divergences) which prohibit the passage from renormalized TMDPDF 
to PDF. This is true whether the soft function is included or not, in the definition of the 
TMDPDF. 

To do so we need the following integrals in d = 2 — 2e and in MS scheme 
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from which we get 
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Taking the sum of the previous equations one gets 
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(4.10) 
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which is the PDF shown in Eq. (4.7). To the best of our knowledge this consistency 
check between the TMDPDF and the PDF is established for the first time. We should 
also mention that even without the inclusion of the soft function in the definition of the 
TMDPDF then one can repeat the above manipulations by considering only the naive 
collinear matrix element fn{h) ^-nd then recover the bare PDF by integrating over transverse 
momentum coordinates. 

Now, from Eqs. (4.1) and (4.3) we extract the following condition for the matching 
coefficient at the intermediate scale in momentum space 

oo 

d'^knL Cn{x; kn±, Q,IJ)= / dfknL X! OilCniix] kn±, Q, n) = 5(1 - x) , 

c^''A;n±C„i(x;A;„x,Q,/^) =0, Vi > , (4.12) 
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where Cnoix;kn±, /J.) = S{1 — x)6^'^\kn±)- It should be stressed that the last equations 
apply to the UV-finite contributions of the integrals of Cni{x; kn±, fi) over kn± since clearly 
those integrals are UV-divergent in momentum space. 

The coefficient at 0{as) in Eq. (4.8), once we transform it back to momentum space, 

is 

Cni{x;kn±,Q,iJ,) = ^ I (27r)^^ 
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\ ^ / ^IR y 



IR eiR fJ' 



Integrating over kn± one can easily see that all the IR poles cancel and the remaining 
contribution contains only single and double UV poles. Actually one gets 



r,/, - ^6{l -x)+ ln^5(l - x)^ ^ - S{1 -x)^ 
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thus the UV-finite term is zero and Eq. (4.12) is established for i = 1. 
5 Q^-Dependence and Resummation 

The matching coefficient Cn is expected to live at the intermediate scale ~ 1 /b. However 
the appearance of ln((5^//i^) in Cni, and higher powers of it in higher orders in perturba- 
tion theory, might indicate otherwise. Notice, for example, that the logarithms in Eq. (4.8) 
cannot be combined into a simple logarithm, unlike the case of threshold region in inclusive 
Drell-Yan or DIS [35, 36]. In the threshold region the matching coefficient at the interme- 
diate scale m is a function of only one logarithm, ln(//j///^). Nonetheless, from general 
arguments concerning the ^-regulator we can extract and exponentiate this Q^-dependence 
in the TMDPDF itself, thus putting it under control to all orders in perturbation theory. 

Working in pure DR and setting all scaleless integrals to zero, only real diagrams 
contribute to jn- Then, we can express the logarithm of the TMDPDF in impact parameter 
space as 



where 



Injn = In/n 
ln/„ = Tin ^x; as, Lt, In 
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and we have set = A. The need for (5-regulator to regulate rapidity divergencies in 
individual Feynman diagrams of fn and 4> introduces the logarithmic dependencies shown 
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in Eq. (5.5). Due to dimensional arguments and Lorentz invariance, those are the only 
possible combinations that can appear. 

Since the PDF is zero in pure DR and the matching coefficient between the TMDPDF 
and the PDF does not depend on the IR regulator, we have 

;^l.i„ = 0, (5.3) 
which implies that TZ^ and TZ^ must be linear in their last arguments. Thus we can write 

lni„ = Inj^' - Dias^Lr) + Lt^ . (5.4) 

where we have introduced Lt just to cancel the /^^-dependence in the coefficient of D which 
simplifies the RG equations of the TMDPDF. The function Inj*"^ is independent of and 
all the Q^-dependence appears explicitly only in the ln((5^//z^). Hence, we can extract all 
the Q^-dependence from the TMDPDF and exponentiate it, putting it under control and 
building what we will call the "Q^-factor" hereafter. 

We believe that the linearity in ln(Q^//Li^) can be extracted without relying on a 
particular scheme of regularization, but based on general arguments concerning the rapidity 
divergencies. As we have shown in Sec. 3 to first order in as, the TMDPDF is free from 
rapidity divergencies, since all the A-dependence that remains exactly matches the IR 
contribution of full QCD. Then, although our (5-regulator does not differentiate the origin 
of the divergencies that it regulates, i.e., it encodes both the IR (soft and collinear) and 
rapidity divergencies, actually one could use another regulator that makes this distinction 
manifest. For instance the z/-regulator introduced in [34]. 

Now, if we denote by v the parameter that regulates only the rapidity divergencies (and 
using a different ones for the IR), then we believe that, based on the 0{as) calculation, 
the functional dependence of ln/„ and liKp on u should be to all orders 

,2 






ln/„ ln-2 



Inc^ — > In— = In— + In— , (5.5) 



where we have taken p'^ = p = Q. 

Since we know that the TMDPDF is free from rapidity divergencies, thus one can write 

:7^1nj„ = 0, (5.6) 
amv 

regardless on how the IR divergencies were regulated. And this equation again implies that 
TZn and TZ^ must be linear in the logs of which automatically leads to Eq. (5.4) and the 
extraction of to all orders in perturbation theory into the Q'^-ia.ciov. 
Using Eq. (4.1) the TMDPDF can be written as 

/g2^2g27B\ -^("-^^) 

jn{x;h^,Q,^) = \ C„(a;;6_L,/^) 8) Qn(a;;^) , (5.7) 
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where 



Cn(x;6_L,/i) = (5(1 - x) + 



27r 



-r,/,LT + (1 - x) - ,5(1 -x){--L't- ^Lt + - 



(5.8) 



The important thing to notice is that all the Q^-dependence in the TMDPDF is exponenti- 
ated to all orders in perturbation theory where the exponent D is perturbatively calculable 
and Cn is Q^-independent. Notice also that Eq. (5.7) refers to one single TMDPDF, and 
not to the product of both as in [8] . 

Given the renormalization group invariance of the hadronic tensor M in impact pa- 
rameter space, 

M = H{Q^/^t^)jn{x;b^,Q,l^i)jn{z-,b^,Q,l^) , (5.9) 

(see also Eq. (7.1) below), we can establish the following relation between the AD of the 
hard matching coefficient, jh, and the one of the TMDPDFs, 7n(fi)) 

IH = -In - 7n = -27n , (5-10) 

where 7„ = 7^ and 

The AD of the hard matching coefficient is linear in ln{Q'^/fj?) to all orders in perturbation 
theory [36, 37], 

-fH = A{as)ln^ + B{as), (5.12) 

where A{as) and B{as) are perturbatively calculable and are known up to third order in 
as- Thus we get 

7„ = --A(a,)ln^--S(a,). (5.13) 

Applying RG invariance to the cross section, and the fact that A[as) = 2rcusp(cKs) to 
all orders in perturbation theory, we get 

dD{as,LT) , . , . 

= rcusp(as) • (5.14) 

am/x 

The perturbative expansion of D is 

D{as,LT) = Y.dn{LT)[^) , (5.15) 



n=l 



where cii(Lr) can be straightforwardly extracted from Eq. (4.6) and it is: di(Lj') = 2CfLt- 
d2{LT) can be read off from the result in [8] by taking half of their result for d\{LT)- The 
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factor of half results from the fact that we are considering only one collinear sector rather 
than a combination of two. Thus 

d2{LT) = ^L't + \^iLt + CfCa - MCs) - (^) CfTfuj , (5.16) 

where we have used the following expansions of the cusp AD and the beta function l3{as) = 
das/dlnfi 

rcusp(a.) = Ern-i(^) , /?(«,) = -2a, ^/3„_if^j . (5.17) 

n=l n=l 

It is worthwhile at this stage to compare our analysis for the dependence of the 
TMDPDF with the one in [8]. There are two major differences: the first is related to the 
origin of this hard-scale dependence and the second is to which quantities it contributes. In 
[8] the emergence of the at the intermediate scale is attributed to the so-called "collinear 
anomaly" resulting from the use of the analytic regularization scheme. This scheme breaks 
the symmetries of the classical SCET Lagrangian and thus it can mediate an interaction 
between two distinct and otherwise decoupled two collinear sectors. In this way the collinear 
anomaly compensates for the absence of a soft function in the factorization theorem for the 
(^T-dependent DY process obtained in [8]. Moreover this collinear anomaly appears only 
in a product of two transverse-momentum-dependent matrix elements. In this product all 
dependence on the analytic regularization parameters (a and /3) is canceled however the 
Q'^ dependence emerges. As a result one can indeed recover the full QCD result for the 
relevant hadronic tensor. In our case the Q^-factor appears for each one of the individual 
TMDPDFs and it is due to the existence of the soft function in the definition of the 
TMDPDF. We also mention that in the analytic regularization scheme the soft function 
vanishes at any arbitrary order in perturbation theory due to scaleless integrals so, trivially, 
only the tree level contribution survives. If one adopts this scheme in our definition of the 
TMDPDF then the remaining naive collinear contribution becomes ill-defined as argued in 
[8]. Turning this question around we might consider the factorization theorem in [8] but 
with a regularization scheme in which the soft function introduced here does not vanish. In 
any such scheme (whether it is the (5-regulator, non-zero offshellness, massive quarks and 
gluons or even going off-the-light-cone) the factorization theorem in [8] will not reproduce 
the full QCD resuh. 

5.1 Resummation 

In the kinematic region where Aqcd Q the logarithms of the scales ratio need 

to be resummed to all orders in perturbation theory. For phenomenological applications 
one also needs to consider the DGLAP evolution of the PDF from a factorization scale 
up to some intermediate scale as illustrated in Fig. (5). The DGLAP evolution is 
well-understood and will not be discussed any further below. In impact parameter space 
where the factorization theorem becomes a simple product one might be tempted, following 
the effective field theory methodology, to resum large logarithms of {qx/Q'^) by evolving 
the relevant anomalous dimension(s) of the effective theory operator (s). This would be 
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Figure 5. Structure of Drell-Yan factorization theorem. QCD is first matched onto SCET-qx at 
the scale ji Q through H , followed by the RG running down to scale ^ ~ qx, resumming part 
of the logarithms ofqx/Q- Then the TMDPDFs jn(n} are matched onto the standard PDFs Qn(n) 
at the scale fi ^ qx, in the impact parameter space, through Cn(n)- The rest of the logarithms of 
qx/Q are resummed by the exponentiation of the -factor. Finally, the PDFs are evolved from 
fi ^ A.QCD up to fi ^ qx via DGLAP equations, resumming logarithms of h.QCD I qT ■ 



true in the case of threshold resummation however this pattern is not sufficient to resum 
ah logarithms for Iow-qt observables. As was pointed out in [8] and as we mentioned 
in the previous subsection, the appearance of the logarithmic Q^-dependent terms in the 
OPE Wilson coefficients-order by order in perturbation theory-of the TMDPDF onto the 
integrated PDF complicates the standard EFT resummation procedure since, on one hand, 
those logarithms do not cancel by any choice of the intermediate scale and, on the other 
hand, they cannot be resummed by standard RGE equations. They are resummed once 
they exponentiate and form the "collinear anomaly" contribution as illustrated in [8]. 

As we have mentioned in the previous subsection, the origin of such logarithms is at- 
tributed, in our case, to the non-vanishing contribution of the soft function to the TMDPDF 
rather than collinear anomaly. However the resummation of the large logarithms can still 
be preformed in the same way as is done in [8]. In both cases we have the hard matching 
coefficient is identical and also the final form of the factorization theorem (after the OPE 
is performed.) It is clear then that the resummation procedure for the hadronic tensor can 
proceed along the same lines. The major difference though, is that in our case it is possi- 
ble to discuss the resummation of large logarithms contributing to individual TMDPDFs 
rather than to the complete hadronic tensor. This fact is important phenomenologically. 
One can obtain a resummed TMDPDF in one high-energy process and implement it in a 
different one due to the universal features of this quantity. More discussion about this is 
given in Sec. 6. 

The resummed hadronic tensor is 



X exp 



Q 



(5.18) 
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where the resummed TMDPDF in impact parameter space is 



jnix;b±,Q,iJ.) = exp 



exp 



^' da' ■ 

—In 
—In 



jnix;b±,Q,ni) 



Cn (x; ///j (gi Qnix; ^j) . 

(5.19) 



All the large logarithms in Eq. (5.18) are contained in the first exponential, the Q^- 
dependent factor and the evolution of the PDF (in Eq. (5.19)). When matching QCD 
onto SCET-gr we extract the coefficient H, and by running it from Q down to /x/ we 
resum part of the logs of qx/Q- The rest is resummed by the exponentiation of the Q^- 
dependent factor, which comes from the OPE of the TMDPDF in SCET-gy onto the PDF 
in SCET-II. Finally, since we need the PDFs at scale /i/, they are evolved by the standard 
DGLAP from a lower scale Aqcd up to fij, resumming all logs of ^qcd/qt- Notice that, 
due to Eq. (5.10), the running of the hard matching coefficient H from Q down to jij with 
7// , is actually equivalent to the evolution of the two TMDPDFs from ^7 up to Q with 7^ 
and 7n. ^ 

The AD of the TMDPDF at first order in was already given in Eq. (3.20). Based 
on Eq. (5.10), we can extract it from [16] at second order in a^: 



7n2 



vr 




18 



-Nt 



+ ('fc(3) 



961 



11 



16 X 27 48 



— 7i' CaCf + 777 + 



65 



24 8 X 27 



NfCp+\- 



^-3C(3)jc|^ 
(5.20) 



The last result and the AD at third order in a^, 7^3, which can be extracted in the 
same manner as 7„,2 from [38] (see also [39, 40]), are essential ingredients to perform phe- 
nomenological predictions with higher logarithmic accuracies. We can write the resummed 
TMDPDF in momentum space as well. 



jn{x] kn±, Q, l-l) 



(2vr)^ 



X exp 



—In 



'g252g27B\-^(°-^^=0) ^ 

J 1 Cni^x;b_i_,fiij^Qnix;ni). 

(5.21) 



Notice that the expression above suffers from the well-known Landau pole when integrating 
over large values of 6, since the integrand depends on as{iii). In the literature this issue is 
generally overcomed by setting a cutoff in b and adding a non-perturbative model function 
for the contribution from long-distance physics. However, when the resummation is done in 
momentum space, following the procedure explained in [8], one would expect to sidestep this 



^Notice that 7„(fi) refers to jn{n) and not to 
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issue for individual TMDPDF. The resummation of large logarithms directly in momentum 
space for individual TMDPDF (and not the product of two) can also be accomplished 
following similar analysis as the one advocated in [8] however this issue will be addressed 
elsewhere. 



6 Universality of the TMDPDF 

The predictive power of perturbative QCD relies on the universality of the non-perturbative 
matrix elements that enter the factorization theorems relevant for different high energy 
processes. Those quantities can be extracted from a limited set of hard reactions and 
then applied to make predictions for other processes. In the following we examine the 
universality properties of the TMDPDFs, Eq. (2.13), by considering them in two different 
kinematical settings: one is for DIS and the other is for DY. 

The difference between DIS and DY settings appears already at the level of the operator 
definitions of the collinear and soft matrix elements of the TMDPDFs due to the existence 
of different Wilson lines between the two settings. Moreover, and since the soft function 
connects two collinear sectors, which are obviously different between DIS and DY, it is not 
immediately clear how the universality of the TMDPDFs is realized. 

For completeness we write below the Wilson lines for DIS kinematics, which can be 
used in actual calculations for the TMDPDF that obeys the DIS setup: 



n{n) n.(n) ^n(ri) i 



Wn{x) = Pexp 
fn{x) = Pexp 
Sn{x) = Pexp 
Tsn{x) = Pexp 



-ig I dsn-An{x + ns) 



-^9 



/ dr l± ■ An±{x^ ,oo~ ,x± + 1±t) 
Jo 



-ig I ds n ■ As{x + ns 

-ig / dr l± ■ As±{0'^ ,oo~ , x± + 1±t) 
Jo 



(6.1) 



where the rest of the Wilson lines can be obtained by exchanging n n and P ■h^ P and 
the relevant matrix elements for the TMDPDF in DIS are 

fn'^={p\ [inW^] (0+,r-,r-l)| \wr^ (0)|p)|,b included, 



f- 

J n 



DIS 



i 



inWl\ (0)^ ^Vin\ (r+,0-,fx)|p)|,b included, 

5Tt oT 



0Z5,5 = (O|Tr Slsl^ (0+,0-,f^) SVSI (0) |0) . 



(6.2) 



Notice that all the collinear Wilson lines are different between DY and DIS, unlike the 
soft ones, where only the n-soft Wilson line changes. This is due to the fact that collinear 
Wilson line, say W^, "knows" about the collinearity of the n-sector, while the soft ones are 
related just to their own sector. 

In this section we will show to first order in ag that the TMDPDF, defined in Eqs. (2.10) 
and Eq. (2.13), is the same in DY and DIS kinematics. This expectation is based on the 
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results of Sec. 4. In that section we have shown that when we integrate the TMDPDF over 
kn± we recover the standard PDF. It is well known that the PDF is universal. Thus, it is 
expected that our TMDPDF will also be universal (at least to first order in as-) Below we 
show this to hold using the 5-regulator and in momentum space. 

The virtual contribution to the TMDPDF in the case of DY is given in Eq. (3.10). 
The real contribution comes from the combination of Eqs. (3.12, 3.13, 3.14) according to 
Eq. (3.6). 

The TMDPDF that enters the factorization theorem for DIS at low qj- spectrum is 



(2^)^ ./0^^^^(O+,O-,f^) 



(6.3) 



where f^^^ and (p^^^ were defined above. The diagrams in figs. (1) and (2) give collinear 
and soft virtual contributions respectively to j^^^ ■ The WFR diagram (la) and its Her- 
mitian conjugate give 



1 1 1 

h -In— ^-r h - 

2euv 2 -iA- 4 



+ h.c. 



(6.4) 



The W Wilson line tadpole diagram, (lb), is identically 0, since n? = 0. Diagram (Ic) and 
its Hermitian conjugate give 



fDIS,{lc) 
Jnl 



-2ig^CF6{l-x)6^^\kn±)ti 



2e 



+ k+ 



27r 



{2tiY [k+ + [{V + + iA-] [k"^ + iO] 
25+A- A- . oA- 



+ h.c. 



2 8+ 2 
-In— H 



In" 



euv P"* 



21n^ + In^^ + 2 H 

12^ 12 



(6.5) 



The contribution of diagrams (2a) and (2b) is zero, since (2a) is proportional to = 
and (2b) to = 0. The diagram (2c) and its Hermitian conjugate give 



d'^k 



{27tY [k+ + i6+][k~ + i6~][k^ + iO] 



+ h.c. 



2-K 



2 , 5+6- , .6+5' 
-In — ^ + ln^ — - 



TT 

T 



(6.6) 



Thus, the results for the virtual contribution for naive collinear, soft, TMDPDF and pure 
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collinear are, respectively, 



fv,DIS _ 
Jnl — 


fDIS,0-c) 
Jnl 




OisCp 


iV,DIS 
<Pl 


j_DIS,(2c) 
91 


■v,DIS 
Jnl — 




fV,DIS 
Jul — 


2-K ^ 



2 Jnl 

<5(l-x)5(2)(4^) 



2 , A 3 1 
-In— + 



euv 2 euv 



In^ 



3, A oA 
2^V + ^ 



7 Svr^ 
^ 4 ^ 12 



27r 



6^^Kkn±) 



-uv 



A2 



euv 



-In 



+ ln^ 



A2 



1 1 

+ 



^uv ^uv \ 2 



2 1 

+ 



^uv ^uv \ 2 



3 , n 
+ 2hi 



— In— :r 
2 /i2 

3, A 



vr 

T 



1. 9 A2 



-In^ 



+ ln^ 



7 

+ T + ^ 



A 



In^ + In" 



A 



+ 



7 tt' 



12 
(6.7) 



Notice that the pure cohinear matrix element given in the last line above is calculated by 
subtracting the complete soft function from the naive collinear matrix element. For DY 
one can get the same result. Comparing our results for DY kinematics with DIS ones we 
get the following 



tv,DIS 
Jnl 



rni + 



2tt 



5{1 - X)d^^\kn±) 



vr 



and 



iv,DIS 



^1 + 



2-K 



(6.8) 



(6.9) 



We show below that the real part of the naive collinear and soft matrix elements are 
also different, and that this difference exactly compensates the one in the virtual parts. 
We remind the reader that all the results below are valid for infinitesimally small A^ with 
respect to all other scales. Diagram (3a) is the same for DY and DIS, 



?{3a) _ fDIS,i3a) _ asCp , w 
Jnl —Jnl ~ ~2tT^ 



x) 



\kn±\'^-iA [l-x] 



2 • 



(6.10) 



The contribution of diagrams (3b+3c) for DY was given before and it can be expressed as 

1 



?(3b+3c) _ asCp 

Jnl - 

~ ~2^ 



(1 - x) +i5+/p+ 



PV 



1 



\k 1 |2 



X 



iA-(l - x) 



+ 



+ h.c. 



{I - x) +i6+ /p+ {l-x)-i5+/p+ 



(6.11) 
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while for DIS it is 



fDIS,{3b+3c) . 2^ 

fni = -"^^9 CfP 



k+ 



X (5 ((1 - - A;+) (5(2) (a;^ + 4^) + /i.e. 



27r2 

asCp 
27r2 



(l-x) -i(5+/p+ 



zA-(l - x)_ 



+ h.c. 



\k 



nl. I 



+ 



{I- x)+i6+ /p+ {I - x) -i5+ /p+ 



+ 



(6.12) 



Thus, the real part of the naive collinear matrix elements in DY and DIS kinematics are 
related by the following 



?r,DIS fir 

J nl Jnl 



f 



nl 



27r2 

2tt 



6{1 - x)6{\kn±\^)27r^ 

d{l-x)5^^\K^)7T\ 



(6.13) 



where we have used: (5(|A;„^|2) = (7r/2)5(2) Combining the last result with Eq. (6.8) 
we conclude that the naive collinear matrix element is universal to 0{as)- 

The soft contribution in diagrams (4b+4c) for DY was given in Eq. (3.14) while for 
DIS we have 

d'^k 



^Dis,{4b+4c) 



(27r)'^ 
1 



5^^\k^ + K±)S{k^)0{k+) 



1 



[k+ - id+][-k~ + id" 



-In 



6+6- 



+ h.c. 



(6.14) 



Since we are interested in expressing our results in terms of distributions in momentum 
space, it turns out to be easier to consider the difference between the real contribution of 
the soft function for DY and DIS. In order to achieve this let us write the following: 



,D75,(46+4c) 



2 9^^{a;t), 



5^^(a;t) = 



1 Int 



at - 1 ' 
1 Int 



at + l 



where t = \kn±\'^/{6+6 ) and a = 6+6 . One can easily see that 

<7^^^(a;t)-5^^(a;t)=A5(2)(4^), 

and integrating over kn± we get the coefficient A, 



(6.15) 



(6.16) 



A 



(fkn± 



DIS 



{a;t)-g^^'{a;t) 



an 



dt 



DY I 



2tt 



~ , Int 
dt- 



TT 

T 



(6.17) 
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The functions and g^^^ are UV-divergent when integrated over kn±. However when 
we take the difference, we get a UV-finite contribution and IR-regularized with the a 
parameter. This is due to the fact that the difference between DY and DIS integrals is just 
the position of the pole of the collinear Wilson line, which is related to the IR (collinear) 
divergence. Thus, the real contributions to the soft functions for DY and DIS kinematics 
are related according to 



■,DIS ir "«^^r(2)/r \ 2 



^1 



2tt 



5^^\knl_)7,\ (6.18) 



Combining this result with Eq. (6.9) we conclude that the soft function is universal to 
0{as). 

To conclude this section, we have shown that the naive collinear and the soft are 
universal, from which the pure collinear and the TMDPDF are clearly universal. In the 
Appendix we calculate the TMDPDF in impact parameter space for DY and DIS, and then 
match it onto the PDF, where all those quantitates are calculated with the 5-regulator. 
By doing so, we show that the PDF is universal, as it should be, and that the matching 
coefficient at the intermediate scale is the same for DY and DIS kinematics and independent 
of the IR regulator. 

7 Factorization at 0{as) 

In this Section we establish the factorization theorem given in Eq. (2.17) to first order in Ug- 
We do it through Eq. (2.14), since we have already established the OPE of the TMDPDFs 
in Eq. (4.1) given the results in Eqs. (4.7,4.8). The hard matching coefficient for the qt- 
dependent DY cross section is the same as the one for inclusive DY. As mentioned before, 
this matching coefficient at the higher scale Q is obtained by matching the full QCD cross 
section onto the imaginary part of the product of two effective theory currents. This echoes 
the "subtraction method" in perturbative QCD. 

We start by rewriting Eq. (2.14) in a more useful way, 

= 3N,q^s 2(2vr)4 «.^(^' ^^^Q) ^ 



M{x,z-q^,Q) = HiQ^fi^) J d^K±d^kn±6^^Hq^ - kn± - 4±) [<5(1 - x)5^^\K±)5{1 - z)6^^Hh 

5(1 - z)6^'^\kn±) + jni (5(1 - x)5^'^\kni_)y\ + 0{al) 
= H{Q^/fi^) [6{1 - x)6{l - z)6^'Hq^) 

+as(s{l - z)jni{x;q±,Q,fi) + 6{1 - x)jni{z;q±,Q,nfj + O(a^) , (7.1) 
where M is the hadronic tensor. 
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The real parts of jn and jn at order as are given in Eqs. (3.15)-(3.16) and using 
Eq. (7.1) one obtains the total real part of the hadronic tensor M in the effective theory, 



^SCET 



27r2 g2 
+5{l - z) 



6{l -x){l-z)+ (5(1 -z){l-x) + 5{l - x) 



2z 



{i-zy 



+ 25{l - x)5{l - z)\n% 

(1 - x)+ 



(7.2) 



The real part contribution for DY at non-vanishing in QCD can be read-off from the 
result given in Ref. [17] (see Eqs. (57)-(59) in that reference) by going from space-like 
(DIS) to time-like (DY) kinematics. Since the QCD result for DIS includes only single log, 
then the results for M*" at 0{as) are equal in both cases. Thus, one can easily see that 



M, 



QCD- 



In QCD the virtual part of M with the 5-regulator is 



M^QCD 



2-K 



6{1 - x)6{l - z)6^^\ 



. 2 A 



A 9 7r2 
'V"2 + Y 



(7.3) 



The above result can be simply obtained by considering the one-loop correction to the 
vertex diagram for qq — t- 7*, with the inclusion of the WFR diagram while using the 
fermion propagators in Eq. (3.1) with = A. 

The virtual parts of j'^ and jn at one-loop are given in Eqs. (3.10, 3.11). Using Eq. (7.1) 
the total virtual part of the hadronic tensor M in the effective theory is 



M^SCET = H{Q^/^?) 



2-K 



5(1 -x)(5(l -z)5(2)( 



-21n^ 



7 



+ 



1 



uv 
27r2 
3 



(7,4) 



where the UV divergences are canceled by the standard renormalization process. We notice 
that the IR contributions in Eqs. (7.3, 7.4) are the same, thus the matching coefficient 
between QCD and the effective theory at scale Q is: 



1 + 



27r 



.1 



8 + 



6 



(7.5) 



The above result was first derived in [35, 41]. We can also obtain the AD of the hard 
matching coefficient at 0{as) and verify Eq. (5.10), 



27r 



6 + 41n 



-27n 



1 • 



(7.6) 



So we conclude that the factorization theorem in Eq. (2.17) is satisfied to first order 
in as- The IR divergences of full QCD are recovered in the effective theory calculation, 
Eq. (7.4). The real contribution is the same in QCD and in the effective theory, and finally, 
the matching coefficient at the higher scale depends only on the hard scale as it should 
be. 
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Figure 6. The soft function at one-loop in light-cone gauge. 

8 Light-Cone Gauge 

In this Section we show that the TMDPDFs, Jn{h)i actuahy the same in Ught-cone 
gauge and Feynman gauge, once the contribution from the transverse Wilson hues is taken 
into account. In Ref. [32] two of us have shown that the naive colhnear contribution to 
the TMDPDF (the numerator in jn(h)) is actually gauge invariant with a one- loop calcu- 
lation. In that article the authors used a particular IR regulator for light-cone divergences 
however the results obtained in covariant gauge and in light-cone gauge are the same and 
independent of that regulator once the zero-bin corrections are included. That was shown 
explicitly in the Appendix of that work. 

In light-cone gauge we use the ML prescription [43], which is the only one consistent 
with the canonical quantization of QCD in this gauge [44]. Moreover in the n and n 
collinear sectors the only gauge fixings compatible with the power counting of the collinear 
particles are respectively nAn = and nAn = 0, which correspond to "killing" the highly 
oscillating component of the gluon field in each sector. We now compare the integrals that 
we have evaluated in Feynman gauge with the corresponding ones in light-cone gauge. 

The interesting contribution to the collinear part of the TMDPDF in Feynman gauge 
is provided by the W Wilson line and it is (cfr. Eq. (3.8)) 



?(lc) {Feyn) 
Jnl 



2e 



p+ + k-^ 



iO) ip + ky + iO ■ 
(8.1) 



In light-cone-gauge this result is reproduced when combining the axial part of the WFR, 



/(^)(^-)=<5(l-x)5(2)(4^)4zg2c7^^2. 



d'^k 1 p+ + k+ r e{k-) e{-k") 



{2t^Y (^^ + ^0) {p + ky + iO 



+ 



k+ + iQ k+ - iO 
(8.2) 



and the contribution of the T Wilson line is 



d'^k 1 p+ + k+ 



(271^ (fc2 + iO) {p + ky + iO 



e{k-) 



1 



1 



k+-iQ k+ + iQ\ 
(8.3) 



It is evident that /^{'^^ i^'^yi^) _ ji^^^) (^) _ ji^") i^'-^) ^he tadpole diagram is null also in 
light cone gauge since the gluon field does not propagate at infinity [32]. 
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In Ref. [33] we have shown that we need the T Wilson hnes also in the soft sector. We 
show this explicitly by considering the virtual corrections to the soft function using the 
gauge fixing nAg = 0. The only one- loop virtual correction in Feynman gauge comes from 
Bg. (2c) 

In light-cone gauge we have two types of contributions: one from the tadpole diagram in 
fig. (6) and the other one is from the T Wilson line. The tadpole contribution in fig. (6b) is 
zero because the transverse gluon fields do not propagate at infinity as mentioned earlier. 
Explicitly, the tadpole contribution from fig. (6a) is 



'1 



-8^-\K^)2^iC,^- J (2v^)^[fc+]M./fc-+.0ifc--^0fc2 + ^0 
-^(^)(4x)2.,^C.^-/^( 



d'^k ( e{k-) ^ 9{-k-)^ 1 k- 1 



{2TTY\k+ + iO k+ - iOJ k~ + iO k~ - iO fc^ + 
= 0, (8.5) 

because when integrating over k~^ all poles lie on the same side of the complex plane. 
Finally the contribution of the T Wilson line in fig. (6c) is 

V nx; y y (27r)'^ Vfc+ - iO A;+ + iO/ fc" + iO F + iO ^ ' 

Notice that we can add to (jii"^'' the quantity /"^ = which is defined as 

/-=-*p)(4x)2«^cp,- /'-^:|j(#cL+ (8.7) 

V n_L; y y (27r)'^ VA;+ + 20 k+ - iO) k" + iQ k^ + iQ ^ ' 



The quantity is exactly zero because when integrating in /c"*" all ^oles, a^ain, lie on 
the same side of the complex plane. Now it is easy to verify that (f)\ '^^ = t/)]^ + at 
the level of integrands. In other words the T Wilson lines in the soft sector insure the 
gauge invariance of the soft matrix element irrespective of any infrared regulator. Similar 
considerations hold for Feynman diagrams with real gluon contributions. 

As a final comment let us consider the work of Ref. [45] where the authors consider 
the definition of TMDPDF in light-cone gauge. In that work it is argued that in order 
to properly define the TMDPDF one should divide the collinear matrix element by the 
soft function (and not by the square root of it.) However in this Section we have shown 
explicitly that the light-cone gauge provides a reshuffling of Feynman diagrams, while the 
final results for matrix elements are the same in all gauges. We have also shown that when 
subtracting the square root of the soft function, as defined in this work, one can achieve 
a separation of UV and light-cone divergences. It is not clear to us how this feature is 
realized in Ref. [45]. 
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9 Equivalence of Soft and Zero-Bin Subtractions 



Let us start this discussion by considering Eq. (3.13) which gives the non-trivial real gluon 
emission to the naive collinear contribution to the TMDPDF. When taking the gluon 
momentum k to the soft limit: k ~ Q{X, X, X), one needs then to distinguish between 
generic values 1 — x where it scales as 1, on one hand, and the threshold region where 1 — x 
scales as A on the other. In the former case, taking the soft or zero-bin limit amounts to 
dropping the k^ from the 5((1 — x)p'^ — k^) thus getting a trivial 5(1 — x) contribution. In 
this case the equivalence of soft and zero-bin subtractions can be easily verified, as we show 
below. However at the threshold region and in the soft limit the term 6{{1 — x)p'^ — k'^) 
remains intact. This will give a non-trivial x-dependence, manifested not only by 5{1 — x) 
but also with the appearance of 1/(1 — x)^ in the zero-bin contribution at 0{as) and 
with more involved "+" distributions at higher orders. Given that our soft function is 
independent of x, then the equivalence of soft and zero-bin subtractions breaks down. This 
is in complete contrast to the case of partonic observables at threshold. In the latter, the 
soft function has to have an explicit x-dependence-which arises from separation of the 
soft Wilson lines in the soft function along one light-cone direction-and this dependence is 
fundamental to establish the equivalence of soft and zero-bin subtractions [23]. 

Moreover, when certain IR regulators are implemented different results for the soft and 
zero-bin contributions are obtained. In [34] the zero-bin is zero beyond tree-level while the 
soft function has non-vanishing contributions to all orders in perturbation theory. Below 
we establish the equivalence of the zero-bin and the soft function subtractions at order Og 
while staying on-the-light-cone and using the 5-regulator. The key point to notice is the 
relation between the regulators in both collinear sectors, Eq. (3.3). 

The pure collinear matrix element fn is calculated by first integrating over all momen- 
tum space and then subtracting the soft limit. Clearly this is done on a diagram- by-diagram 
basis and perturbatively, 



fn = /no + (/nl " fnl,zb) + 0(a?) • (9.1) 



We show below to 0{as) that this can be achieved by dividing the naive collinear matrix 
element by the soft function 

fn = J = fnO+ (/nl " fnO <^l) + 0{al) . (9.2) 

For fn analogous analysis trivially applies. The zero-bin of the WFR in diagram (la) is 
zero, and also the one for diagram (lb). The zero-bin of diagram (Ic) is obtained by setting 
the loop momentum to be soft, k ~ (A, A, A) 

= -6{1 - x)5^^\K^)2ig^CFfi^' f ^-^TTT , T + h.c. 

^ ^ V n-Ly y i-p y {2TT)'i[k+ -i6+][p+k- +iA-][k^ + iO] 

= 6(1 - x)0p'^ . (9.3) 

It is clear that with the relations in Eq. (3.3), the subtraction of this zero-bin is equivalent 
to dividing by the soft function in Eq. (3.9), proving the equivalence at order as for the 
virtual contributions. 
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Let us now consider the real diagrams in fig. (3). The zero-bin of the diagram (3a) is 
zero. Diagram (3d) and its zero-bin are zero due to = 0. The zero-bin of diagram (3b) 
and its Hermitian conjugate (3c) is 

+ h.c. 

= 6(1 - x)0S^^+4c) (3 4) 

which is equivalent to divide by the soft function diagram (4b) and its Hermitian conju- 
gate (4c), given in Eq. (3.14), thanks again to the relation in Eq. (3.3). In conclusion, we 
have proved that subtracting the zero-bin is equivalent to divide the naive collinear matrix 
element by the soft function to first order in Og. 



10 Conclusions 

In this work we have studied the Drell-Yan lepton pair production at moderately small 
transverse momentum qt- The analysis was carried out through the framework of the 
effective field theory via successive a two-step matching procedure: QCD — >■ SCET-g^ 
SCET-II. We established an all-order factorization theorem which allows for a phenomeno- 
logical study of DY qt spectrum to be analyzed at energies much larger than Aqcd- When 
considering the double-counting issue of the soft and the naive collinear regions properly, 
the obtained factorization theorem serves as a guideline towards how the TMDPDF should 
be defined and what would be its fundamental properties. In our calculations we have used 
Wilson lines defined on-the-light-cone and light-cone singularities appearing in individual 
Feynman diagrams are regularized with the (5-regulator. We have also introduced our result 
for TMDPDF using pure DR. Based on the relations in Eqs. (2.10, 2.13, 2.14), we were 
able to define an on-the-light-cone TMDPDF which has the following novel features: 

• It can be integrated over the transverse momentum (of the parton in a parton) to 
recover the partonic PDF. 

• It is free from rapidity divergences/logarithms. 

• Its evolution is governed by a single parameter RG evolution equation. 

• It is defined in a gauge invariant way among regular and singular gauges. This 
definition is obtained from first principles of SCET. 

• It is universal among the DIS and DY kinematics and it can be readily modified to 
obtain its gluonic version relevant for the proton-proton high-energy collisions. 

The inclusion of the square root of the soft function in the definition of the TMDPDF 
has an important consequences: 

• The double counting in the factorization theorem is taken into account. 
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• It allows for the separation of UV and IR divergences in the TMDPDF. 

• Even with the subtraction of the square root of the soft function we are able to recover 
the PDF from TMDPDF. 

Staying on-the-light-cone has a set of advantages over going off-the-light-cone: 

• The evolution of the TMDPDF is governed with a standard /_f-RGE. 

• All the results are well-defined without need of using additional parameters. 

• It is consistent with the power counting of SCET. When going off-the-light-cone one 
needs to introduce a tilted vector v{v) where there is no specified power counting (in 
terms of A) among its large and small components. 

The two step factorization is necessary to perform the resummation of logs of qt/Q 
and A.qcd/qt respectively and we have discussed the resummation procedure in impact 
parameter. We also commented on the resummed TMDPDF in momentum space. In 
the first step of the factorization one gets the usual structure of the cross section given 
in Eq. (1-1). The matrix elements so defined however are not good objects for the 
second factorization because of the presence of mixed UV/IR and rapidity divergences. All 
these divergences however disappear in the TMDPDF as we have defined it. The absence 
of rapidity divergences allows the resummation of all logs without the Collins-Soper like 
evolution equations. The second factorization is built up by matching the TMDPDF onto 
the PDF for large qt- We have studied the Wilson coefficients that appear in the second 
matching which contain, in impact parameters space, ln(Q^/^^). We have shown that this 
kind of logs can be exponentiated in a similar way as was done in [8] , although in our case 
the concept of collinear anomaly is irrelevant due to the existence of the soft function. 

We also considered the TMDPDF with DIS kinematics and pointed out the differences 
with respect to the DY ones. As mentioned earlier, different Wilson lines are needed for the 
two kinematical settings. However we established the universality of the TMDPDF in both 
regimes and argued its validity to all orders in perturbation theory. The gauge invariance of 
the TMDPDF was also established by computing it in light-cone gauge with the inclusion 
of transverse Wilson lines. The fact that all Wilson lines are defined on-the-light-cone only 
facilitates this computation. We finally comment that our results can be extended to the 
Higgs boson production at low qt as well as to spin-dependent non-perturbative hadronic 
functions, like Sivers and Boer-Mulders ones. The extraction of the resummed TMDPDF 
from HERA and LHC data, as well as its use to make phenomenological predictions is left 
for a future consideration. 
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A Matching the TMDPDF onto PDF with 5-Regulator 

In this appendix we calculate the matching coefficient of the TMDPDF onto the PDF at 
the intermediate scale using the (5-regulator, and show that, as expected, it does not depend 
on the particular choice of the IR regulator and we get the same result as with pure DR. 
We also show that this matching coefficient is the same for DY and DIS kinematics, thus 
establishing its universality to first order in ag, since the PDF is universal. 

A.l Drell-Yan Kinematics 

The virtual diagrams for the PDF are the same as for the naive collinear matrix element 
that enters into the definition of the TMDPDF, Fig. (1). From Eqs. (3.7, 3.8) we get 



(la) ^ ^f^Lsil 



27r 



euv A 



1 

+ 2 



(A.l) 



and 
2(1^=) 



2-n 



5{l-x) 



2 A+ 2 
-m^^ H m 



, A+ A- 
21n— rln^ 



A- Ttt- 
21n^ + 2 



12 
(A.2) 



The real diagrams are the same as in Fig. (3), from which we get 



Q(3«) = 27Tg'CFp^ 



5{k^)e{k+) 
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(1-x) 
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5 ((1 - x)p+ - k 
(A.3) 



and 
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(A.4) 
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where we have used MS-scheme, /U^ 
5+/P+ < 1, 



/i^e'^^/(47r), and the following relations when 
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_ ^ 2x{l-x) ^ 2x _ _ ^, ^ 
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Combining the virtual and real contributions we get the PDF to first order in 



(A.5) 
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(A.6) 



The virtual part of the TMDPDF in momentum space was given in Eq. (3.10), and in 
impact parameter space it reads 



2tt 

-—In — ^ 
2 /x2 



1 1 

— + 

uv 



+ In 
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(A.7) 



The Fourier transform of diagrams given in Eqs. (3.12, 3.13, 3.14), while keeping the A's 
to regulate the IR divergences, are 



|(3a) asCp,^ 4e-27^^ 
fni = ^^(l-^)ln- 



27r 



A-(l -x)62 ' 



(A.8) 



?(3b+3c) _ asCp 
•^"1 ~ 27r 
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TT 



24 2 Q2 



and 



In the above we have used the following identities in d = 2 — 2e: 
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and also 

^ ^ 1 4p-27B 

ifcx-bx^ = vrln— , 

|A:^|2-iA2 -zA262' 

i"'^' fcFTII = ""W• 
when A — )■ 0. 

Finally, setting = A, the TMDPDF in impact parameter space to first order in Ug 

is 

in = Qn + ^ [ - LtV,,^ + (1 - x) - 5(1 - x) (^-Ll - hr + Lrln^ + 

'(A.13) 

where Qn is the PDF given in Eq. (A. 6) and the remaining part exactly equals the OPE 
matching coefficient calculated in pure DR given in Eq. (4.8). 

A. 2 DIS Kinematics 

For DIS kinematics the operator definition of the PDF changes, as we showed in Sec. 6, 



zb included 

(A.14) 



where x = ^nCn and is the collinear Wilson line defined in Sec. 6. In the following 
we show to first order in that the PDF is universal, as expected, although its operator 
definition changes for DY and DIS kinematics. 
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The real diagrams are the same as in Fig. (3), from which we get 
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Combining the virtual and real contributions we get the PDF to first order in 
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which is the same as in Eq. (A. 6) for DY kinematics. 

The virtual part of the TMDPDF for DIS kinematics is 
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The real diagrams and their Fourier transforms are 
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In the above we have used the following identity: 
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when A —7- 0. 

Finally, setting A^ = A, the total TMDPDF in impact parameter space to first order 
in as is 
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where Qn is the PDF given in Eq. (A. 6) and the remaining part exactly equals the OPE 
matching coefficient calculated in DY kinematics. 



B Quark Form Factor in Full QCD with (5-Regulator 

In this appendix we calculate the Quark Form Factor (QFF) in full QCD with the 6- 
regulator for DIS and DY kinematics. We will show that the real part of the QFF is the 
same in both regimes, and that the imaginary part changes sign. We set A^ = A. 
For DIS setting the vertex correction is 
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where we have taken the large Hmit and only the terms proportional to 7". The WFR 
is the following, 

^ ^ i {2nY [{p - kY + iA][k^ + id] ^ 4tt \e 2 -iAj ^ ' 

The final result for the QFF calculated in QCD, with the 5-regulator and for DIS kinematics 
is 

< Pn\r\vn = v-'^^' - Y-^il - V-^il = 

I if I If 
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The analogous results for DY setting are the following, 

y J (2^)d [(^_^)2+^^-][(^+^)2+^^+][^2+^0] 

and the same WFR as in DIS. The final result for the QFF is then 
<Pn\J \P^> =V ^ - 2-^7x - 2 = 



71 



2 



1 H -21n^— - 31n - 2mln— \ 

47r I Q2 _Q2 Q2 2 2 



(B.5) 



In conclusion, the real part of the current calculated with the ^-regulator is the same 
for DY and DIS kinematics. This is due to the fact that our regulator, iA, is imaginary, 
so changing from space-like (DIS) to time-like (DY) does not change the real part. 

Notice that Eq. (B.5) is consistent with Eq. (7.3) when adding its Hermitian conjugate. 
Combining Eq. (B.3) with its Hermitian conjugate and the virtual contributions Jn^^^ 
and jf^^^^ in Eq. (6.7) we recover the hard matching coefficient for DIS between QCD and 
SCET (Eq. (20) in [36] plus its Hermitian conjugate). 
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